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Abstract 


The aim of the present work is to develop a computer program to compute the 
performance parameters like torque, power and efficiency of a Pelton Wheel with 
a given bucket surface configuration. The water jet striking the Pelton Wheel 
buckets is modeled as a series of particles. The force and the torque exerted by 
every particle, both on impact and subsequently during the motion on bucket, are 
computed, using Newtonian dynamics, at every instant of time and integrated over 
entire time span until they lose contact from the bucket. The time history of all 
particles is then used to obtain the total torque obtained by the water jet, as a 
function of time. The bucket shape is represented in local spherical coordinates 
r = f(9, <f>) for which data are supplied on a coarse empirically measured xyz co- 
ordinate grid system. The performance parameters for two Pelton Wheel buckets 
are calculated with the present model and their results are compared. The aim is 
to help in the optimum design of Pelton Wheel buckets, both in terms of operating 
conditions and bucket surface profile 



Chapter 1 


Introduction 


The purpose of a hydraulic turbine is to convert kinetic and potential energy of 
water into mechanical work. On the basis of the available head, three principal 
types of turbines - impulse turbine, the radial or mixed-flow reaction turbine, and 
the axial flow, or propeller, reaction turbine can be used. 

Among the impulse water turbines, the Pelton turbine also called the free jet 
turbine, is the only type being mostly used these days. It is a high head and 
low flow turbine. The head of the supply water is converted into one or more 
high-velocity free jets which are directed against buckets mounted on the rim of 
a wheel. The impact of water on the surface of the bucket and further motion on 
the surface of the bucket produces a force which causes the wheel to rotate, thus 
supplying the mechanical power to the shaft. 

In general, each bucket of Pelton Wheel consists of two nearly hemispherical 
cups separated by a sharp edge at the center. The water jet impinges at the center 
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of the bucket and is divided by the sharp edge, without shock, into two parts 
moving sideways in opposite directions. The water after impinging on the buckets 
is deflected through an angle of about 165° (instead of 180°), so that it may not 
strike the back of the other buckets and retard the motion of the wheel. The change 
in the momentum due to this change of directions transmits an ’impulse’ torque 
to the bucket which translates to work output. 

The efficiency of the wheel is given as 

work done on the bucket 

Tj = 

energy supplied to the bucket 

The efficiency of the Pelton wheel is generally high (>90%). Therefore if 
improvement has to be made, we need to have analytical methods that can detect 
the effects of subtle changes in the bucket shape. Since the shape of the bucket 
plays a very important role in creating the torque, if the bucket surface is not 
proper the kinetic energy and momentum of the jet is not transferred efficiently to 
the bucket. An improper surface may cause too great a loss of kinetic energy as 
the jet impacts the bucket, and/or allow too much energy to leave with the water 
exiting the bucket. So the optimum surface and the optimum operating condition, 
a combination of which gives maximum efficiency, is required. In most text books, 
two dimensional analysis is presented using velocity triangles for the calculation 
of work done on the bucket. This simple analysis is not adequate since it does not 
consider the surface geometry, except the deflected angle. Through experiments 
we can also analyze the effect of the bucket surface on the power developed by 
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the wheel. But changing bucket parameters and doing experiment requires a lot 
of time as well as resources. 

It would therefore, be of great help to the designer if a computational method 
can be found to theoretically evaluate the efficiency and other performance pa- 
rameters of a proposed bucket shape. The speed and relatively low cost of compu- 
tational solutions will allow the designer to evaluate many possible bucket shapes 
before actually constructing prototypes for experimental testing. This problem, 
though an important one, has virtually no published work. The commercial im- 
portance of the problem has ensured that any research on this problem is kept 
secret. 

The general problem can be stated deceptively simply: Given a particular 
bucket shape, can its performance in terms of efficiency, torque, force and power 
generation be predicted? The problem is a very difficult one. Any attempt to 
obtain answers through the traditional routes of fluid mechanics seems to be cur- 
rently impossible- for the jet/bucket interaction involves a highly unsteady, even 
turbulent, flow in a complex geometry involving free surfaces. To attempt to solve 
for flow using Euler or Navier Stokes equations seems to be completely outside 
the scope of present-day computational fluid dynamics. 

In this thesis a method, far simpler and probably more likely to succeed, is at- 
tempted. The jet is treated as a steady stream of perfectly plastic (i.e, deformable) 
“particles” which strike the bucket and slide on its surface imparting momentum 
and energy to the bucket during their impact and subsequent motion on bucket sur- 
face. The particles are assumed to be non-interacting (i.e the motion of each such 
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“particle” is taken to be independent of others) and the motion of each particle 
is computed using the Newtonian laws governing particle motion. The equations 
governing the particle motion are coupled second-order ordinary differential equa- 
tions (much more numerically tractable than the partial differential equations of 
fluid flow) which can be solved by the Runge-Kutta method. However, consider- 
able difficulty is introduced into the problem due to the generally complex bucket 
geometry involved. Once the time history of each particle is known from the so- 
lution of the dynamical equations, the effects of all particles can be integrated 
to obtain the total picture regarding the efficiency, torque etc., generated by the 
bucket. 

In the present work, we have followed the work of Ajay Kumar (1999) to de- 
velop a set of computer programs to analyze the performance of real-life working 
Pelton Wheel buckets. If the number of buckets in the proposed design are not 
fixed then the simulations are done for single bucket, but the program can simu- 
late more than one bucket sharing the jet at a time. The torque imparted to the 
bucket is computed at every instant of time and integrated over the entire time 
span. The geometry of the bucket is usually very complex and is very important 
in the calculation of torque. Therefore the surface of the bucket is represented 
by a fine grid in spherical coordinates (r, 6, <j>) having uniform intervals of 6 and 
4, with the value of r is specified at every grid point. Bi-linear interpolation is 
used to find the interpolated value at any intermediate point. We can change the 
parameters of the bucket and see the effects on the torque developed and study 
of the other parameters to tell which part (configuration of system or the bucket 
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surface profile or both) is to be improved/optimized. One test case of 3 D planner 
bucket is simulated with this program for the validation of the code. 

This work was partly supported by a project sponsored by BHEL Bhopal. 
Performance parameters for two Real life working buckets with surface profiles 
and operating conditions are provided by BHEL Bhopal are also calculated. The 
effect of changing the various input parameters on the performance of tire Pelton 
Wheel is studied with the program. The program basically has two parts, one 
for grid generation, which represents the surface of the bucket in the required 
coordinate system from empirically measured data, and other solves the governing 
equation over the grid representing the bucket surface to obtain the performance 
parameters. 

This work is the culmination of the work which was started by Ajay Kumar 
(1999) and Gaurav Dar (2000). However, only here have all the necessary ele- 
ments of the formulations developed by these previous authors been put together 
to compute the efficiencies etc. of real-life Pelton Wheel buckets provided by 
BHEL Bhopal. 
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Chapter 2 


Theoretical Formulation 


2.1 Introduction 

As already discussed, because of inherent difficulties in the fluid mechanics ap- 
proach to the detailed analyses of flow in a Pelton Wheel buckets, a particle dy- 
namics approach is adopted here. In this the water jet striking die Pelton Wheel 
bucket is modeled as a series of particles. The momentum and energy transfer 
from these particles to the bucket is then computed by using the principles of 
momentum conservation and of Newton’s second law applied to the individual 
particles. In this chapter, the kinematics and dynamics of the fluid particle in 
three dimensional (3D) situations is formulated. First we consider a 3D model for 
calculating the force and torque imparted on the rotating bucket by a single parti- 
cle. A general shape of the bucket whose equation is given in polar coordinates is 
considered. Later, this model is extended for a series of continuous particles and 
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the concept of mass flow rate is introduced. 


2.2 Kinematics and dynamics of a single particle for 3-dimensional 
bucket 

Following Ajay Kumar (1999), let us consider a bucket of general shape rotating 
with a constant angular velocity ui about the A’-axis as shown in figure 2.1(a). 

There are two frames of reference considered, one fixed frame XYZ centered at 
O, while other one moving frame xyz centered at o'. The frame xyz is attached 
to an extended point (o') of the arm whose distance from the origin of the fixed 
frame (i.e from O) is R. In figure 2.1(b), the bucket is represented in spherical 
coordinates which are fixed to the origin o' of the xyz coordinate system. The 
surface equation of the bucket is assumed to be ip(r, 9,(p) = 0 and its gradient is 
given by 


_ 7 * dip „ Id'ip % 1 dip 

Vt/> = — —r q — 9 q 

V dr r 89 rsmOdtp 


8 


( 2 . 1 ) 


Initially at t = 0, the arm is assumed to make an angle % 0 from the Z-axis, a 
particle p b traveling along the water jet parallel to the Y axis strikes the bucket at 
point p b0 with velocity Vj. The position vector of p b 0 with respect to the moving 
(r, 9 , <p) frame is rr + 99 + <p<p . After striking the bucket at point p b0 , the particle 
is assumed to travel along the bucket. At that instant, absolute acceleration of 
the particle is along —V, which is normal to the bucket surface (i.e. the bucket 
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water jet 


vane 




(b) Bucket in spherical co-ordinates 


Figure 2.1: Schematic diagram of bucket in 3D 
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is assumed to be frictionless so the particle is subjected only to a normal force 
and experiences a normal acceleration). According to the laws of dynamics, the 
absolute acceleration of the particle p b is given by 


a abs = a 0 + wx(u)xf)+2wx (r r 4- r88 + r sin 8epj>) + a r f + a&8 + a^ep (2.2) 

where a 0 is the acceleration of point o’ 
a 0 = —ui 2 Rk 

and dr, ag, a $ are the acceleration components of the point p b0 in f, 8, <P 
directions, respectively. These components are given as 
a T = f — rd 2 — rep 2 sin 2 8 
ag = r 'O + 2 fO — rep 2 sin 8 cos 8 
a, $ = rep sin 8 + 2 rep sin 8 + 2 rdep cos 8 

The relation between the unit vectors in the Cartesian coordinates and spheri- 
cal coordinates for the moving frame are given as follows 
f — sin 8 cos epi + sin 8 sin <f>j + cos k 
8 = cos 8 cos epi + cos 8 sin epj — sin 8k 
ep = — sin epi + cos epj 
The above relations can also be written as 
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f \ 

i 


sin 8 cos <p 

cos 6 cos ep 

— sin <p 


f \ 

f 

3 

> = 

sin 9 sin ep 

cos 6 sin p 

cosep 

< 

8 > 

k 

K J 


cos 9 

— sin 9 

0 


P 

\ / 


Also the relations between the Cartesian coordinates of the fixed frame and 
the moving frame are as follows; 

/ = i 

J = cos xj + sin xk 
K — — sin xj + cos xk 
equation (2.2) can be written as 


Clabs — A r f + Agd + A^Cp 


(2.4) 


where 


A r = — uj 2 R cos 6 — w 2 r(cos 2 0cos 2 (p + sin 2 <p) — 2wr(<j>cos#sin 9 

cos <p + 9 sin <j>) + (f — rO 2 — rep 2 sin 2 8), (2.5) 

A$ = uj 2 R sin 9 + w 2 r sin # cos # cos 2 ep + cur sin ep + 2us (rep sin 2 8 

+r sin <jj) 4- (r8 -f 2r6 - rep 2 sin 6 cos 9 ), (2.6) 
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and 


,4<p = — oj 2 r sin 8 cos (psin <j> + wrcos^cos <p + 2uj(r cos 9 cos <p — 

— rO sin 0 cos <j>) + (rcj) sin 6 + 2f0sin# + 2r(9pcos$). (2.7) 

However, as the bucket is assumed frictionless, the net force on the particle, 
and hence its absolute acceleration, can only normal to the bucket surface, i.e., 
in the direction of — Vt/j.This implies a relationship between the magnitudes of 
A r , A e , A,p such that die resultant acceleration is in the normal direction, following 
two equations are thus obtained. 


Ar _ Ae 

1 dv 
dr r 06 


and 


A r _ A$ 

dv 1 dv 

dr r sin 6 dc p 

along with 


( 2 . 8 ) 


(2.9) 


i>(r,e,<f>)= 0 ( 2 . 10 ) 

Substituting equation 2. 10 in equations 2.8 and 2.9 and using expressions (2.5, 
2.6 and 2.7) we get 
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( 2 . 11 ) 


G 110 4 Ol20 4 &13V 4" &l4<p 2 4 Gl50 + CL ie#0 4 c — 0 


and 


021$ 4- 0220 4- 023$ 2 4 O240 2 4 &25^ 4 O26$0 + d — 0 
where all a’s are function of geometric variables and are as follows 



®12 = 


dip dip 
1 d $ 



Ol3 



d l y >dy \3 
’ ' BQ ' 


i Ml 

r (¥ 3 


d z ip dtl? 
1 flflS ~M 


¥ 


_L 

■*" dd 


a 14 


d z v dtp By 
o 1 oroo ~m ~oc 

~~ — — 


—¥(¥* 


& 


4rsin#cos#|^ 


9 i 

jL dr 
r (&) 
)d\h 


d z v dip 

1 

r~W~ 


sin 2 ^f 


ai 5 = —2a; cos 0 sin 0 cos - 2a;r sin 2 0 cos (j)^ + 2a; sin (j> 


dr 


d<p 


d^ip dip dip d'^ v i dip \ 2 d^ip dip d*ip < dtp \2 dip 

2 ar£? "d$ dp i 2 5 rc)V^ dti ' 2 dvdc ~dW 2 Br‘^ ' dO ' '&<j> 

u 16 7 ( Bv Y2 * T ( dv \ 2 7 57 r 77v7 ~ 

' "<Ur" ; 


r m 2 

! O ^ 

' r2 00 


3? 


(^) 3 


c = — ^|^(w 2 i?cos#4a; 2 r(cos 2 #cos 2 04sin 2 0)) — |^(u; 2 .Rsin$ 


(2.12) 
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+w 2 r sin 9 cos 9 cos 2 p -f wr sin p) 


Q 21 = 


6-0 dv 

1 6# «7P 

r sin 0 

dr 


a 22 ~ 


(H) 2 


r sin 0 


(70 ' 


rsinflf 


a 2 v op ao 

n 2 o^' r cTo 

«23 — rsin0 (|£y2 



d z i> dil) 

1 jggPS 
rsin0 

or 


1 00 
sin0 


fl24 — 


o v / dy \ 2 
. drop ^ op ^ 


0^0 tdv ^3 


’sin0 (£E)2 rsin0 (^-) 3 rsin0 


a 2 0 dii> 

^^ + sin0fi 


a 25 = + 2w cos 0 cos p|| + 2wr sin # cos p|^ 


a 26 


2 »H) 2 


r sin 0 

a 2., 


(|*) 2 


+ 


#-* 


rsin 0 


cry d-d (ddy 2 

dr W’\’Sp) __ 0, r pnQ #00 

TcTpT^ y $ r 


ap 


dv 
1 0~P 


( ^) 3 


+ 2 sin 0^ 


a~a oo 

d&dp 0p 


r sin 0 


*57 


2 

rsin 0 


1 00 


(— co 2 Rcos 9 — w 2 r(cos 2 6* cos 2 p + sin 2 p)) — ^(—co 2 r 


rsin0 d<p 

sin 0 cos p sin p + cur cos # cos p) 


Equations (2. 1 1 and 2.12) governing the motion of the particle on the bucket 
surface, are two coupled second order non-linear ordinary differential equations 
in 9 and p. These can be solved simultaneously as initial value problem, and will 
require four initial conditions in terms of 9 0 , p 0 , 9 Q and p 0 (i.e. the values of 9, 
9, o and p at t = 0). The value of 9 0 and p 0 can be calculated by detennining 
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the point at which the particle path along the water jet meets the bucket surface. 
The methodology to obtain 0 O and <p 0 will be discussed in detail later on. Once 
the particle’s initial position (9q,<Po) has been obtained, then to find the values of 0 
and 0 at t = 0, we proceed as follows: 

As discussed in the previous section, after the initial impact the component 
of jet velocity Vj along Vyj is same as the velocity of the point p h0 on the bucket 
surface in that direction but the velocity of the jet is unchanged in the direction 
normal to Vip, i.e. tangential to the surface. This is because no force could be 
applied tangential to the surface of the buckets (which is frictionless), thus the 
tangential momentum of the particle will remain unchanged. 

Before impact, the velocity of the jet tangential to the surface is 


(v J -.V'0)V'0 

|Vt/f 


( 2 . 13 ) 


and the tangential velocity of the particle after impact is 

(v p .V^)V^ 

Vp IWf 


( 2 . 14 ) 


where v p is the total velocity of the particle after impact, 

As the tangential velocity of the particle remains unchanged after impact, 


v i 


(vj.Vip)V'ip 

|V ^| 2 


= v p - 


(v p .V^)Vt/) 

|V^1 2 


( 2 . 15 ) 


and tire normal velocity of tire particle after impact is equal to tire velocity of tire 
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bucket at the point of impact, i.e., 


Vp.V-0 = v pfe .Vt/> (2.16) 

where v pi , = u> x (Rk -f rf) is the velocity of the bucket surface at the point of 
impact. We can also express v p as : 

v p = u) x (Rk + rf) + rr + rdd + rsinOqxp (2.17) 

where the first expression on the right is v p b and the remaining terms are the rel- 
ative velocities (wrt point pb) in the r, 6 , tp directions. Substituting the above ex- 
pression for v p into equation (2. 15) and using equation (2. 16) we get 


(vj.V^)Vt/) 

v iWF~ 




= u>x (Rk+rf) + rf+r96-\-rsin6(f)(j) (2.18) 


Now 


Vj — vjcosxj + VjSinxk 


(2.19) 


Replacing the unit vectors i,j,k by r , § : 4> in equation (2. 19) by using (2.3), and then sub- 
stituting into equation (2. 1 8) and equating the two sides component-wise we get. 
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9 o = 


and 


Rcosdsincp - jrff + + ■v j (cosxcos6sin<j) - sinxsinO) 


jjrsino 


( 2 . 20 ) 




—ioRcoscj) - + Vjcosxcoso - urcosOcosQ 

rsinO 


( 2 . 21 ) 


In this manner we arrive at two of the four initial conditions required for tire 
particle path computation. 


2.3 Particle Force and torque computation 

Equations (2.11) and (2.12) which govern the motion of each particle on the 
bucket surface are two coupled second-order non-linear ordinary differential equa- 
tions, and are solved as initial value problem which require four initial condition in 
terms of 9, 6, (f>, <j>. All the four quantities (6, 9, <j>, d>) at the subsequent time steps 
(starting from the impact time t pi to the time the particle leaves the bucket) are cal- 
culated by solving the equations (2.11 and 2. 12) using a fourth order Runge-Kutta 
(RK4) solution method, the change in 9, 9, <j), $ with time is- used to calculate 
the velocity and the acceleration at every time step. The particle is considered 
to have left the bucket whenever the computed a n (normal acceleration) becomes 
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negative or the particle reaches the edge of the bucket. Using these values, we 
get tiabs (absolute acceleration) from the equation (2.4). The force exerted by the 
single particle is given by 


F = -m p a abs (2.22) 

where, m p (= rh p At ) is the mass of the particle, fh p is the mass flow rate. 

The initial impulse exerted by the particle during the first impact at p h0 is 
obtained as 


n m?[ iwi |Vtf| j 


and the average force is taken to be 


m p ( 


Vj.Vip %,.V'0 

iwi m\ 


) 


(2.23) 


(2.24) 


where v pb = to x (Rk + rf) 

Momentum is imparted by the particle when it strikes the bucket and also dur- 
ing the subsequent motion until the particle loses contact with the bucket. Torque 
imparted about X axis due to impact by a single particle having mass m p . 

M initial = ®Pb X (2.25) 

where opt, = Rk - f rf = ( R cos 9 + r)r - R sin 69 
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Torque imparted by the particle about the X axis during the motion on bucket 

M At = rn p (dpbXa p ) = m p (-Rsin9A ( t,f-(Rcos9+r)A < p9+(A g (Rcos8+r)+A T RAn 8)q>) 

(2.26) 

The total torque due to motion of particles on the bucket is 

M At = m P Y^{oPb x a p ) (2.27) 

So die net torque= X Ma t + X) M initial, where the summation is over all 
particles present in the bucket. 

2.4 Initial positions of particles 

Since particles strike one by one and keep on sliding along the bucket till they lose 
contact from it, the computation for each particle is done separately from different 
initial conditions. For a particle on the center line of the jet, the initial position 
where the particle strikes is: 

X c = 0 

Y c = R sin x + z sin % + V cos X 
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Z c — H = R cos x + z cos x — V sin % 


where, (X C: Y c , Z c ) is the point of impact with respect to the fixed fran ie ’ 


and 

and 


X=(Xo + wf ) is the angle shown in figure (2.1) which varies from xo t0 

P in th e 

(, x c , y C i z c ) are the coordinates of the point of impact on the bucket surfe 06 

• hove the 

moving coordinates system (figure 2. 1 ), and H is the height of the jet w u 
pivot center O. 

It is however unlikely that particle will come only along the center 111 

■ i-ne alnt^S 

the jet has some finite radius l 0 . So we assume that the particle may com & 


a point randomly located on the circular cross-section of the jet. We 
assigning a random position to the point by the prescription: 


do tw s W 



Figure 2.2: Water jet 


1 = lo\/rarH 
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3 = 2nran2 

where rani and ran2 are two different uniformly distributed random numbers 
having values between [0,1] , the random point at any l and j3 are shown in hg ure 
( 2.2 ). This prescription is used to ensure that the points are randomly distributed 
on the cross section. The initial position of the particles will be perturbed from 
tire center line by : 

AX = / cos/3 
A Z = l sin 3 

where A A” and A Z are the variation in the X-direction and the Z-direction, re- 
spectively, from the center line. So now the equations for the initial positions 
are 


X p = AX ( 2 - 2g ) 

Y p = R sin x + sin % + y cos x (229) 

Z p = H + A Z = R cos x + z cos x ~ U si n A (2.30) 

where (x. y, z) is the initial point in the moving frame which has to be determined 
for a given A A', A Z, x, H and R. In these equations, X p and Z p are known pnon, 
and Y p has to be determined in terms of the unknown coordinates y and 2 which 
should also satisfy the equation of the bucket surface (chapter 3) 
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lK r > 9, 4>) = 0 


(2.31) 


Thus equation(2.31) has to be converted to a form Y = f(X , Z) and solved to 
obtain Y p . This requires a root finding algorithm which is explained below. 

The equation of surface tp(r, 9, (j>) = 0 can be cast as 

Tfi(x,y,z) = 0 (2.32) 


Using equations (2.28), (2.29) and (2.30) we can rewrite the equation ^(x, y, z ) 
0 as 

Tp(X p ,Y p ,Z p ) =0 (2.33) 

In equations (2.33) X p and Z p are known, and Y p is the root of the equation. It 
can be solved using Bisection method. In equations (2.28 to 2.30) we first take a 
guess value of y, then z is calculated with this guessed value of y using equation 
(2.30). After getting rr, y, z the corresponding r, 9 , q are calculated as 
r = yfx 1 + y 2 + ~z 1 
6 = arccos (z/r) 

(ft = arctan (y/x) 

The bucket surface is represented as t/>(r, 9, cj)) ~ r - f(9 : 4>) — 0 and the 
data are supplied in the form of a rectangular grid (M x N) in 9 and 6 at each 
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grid point of which the radial distance (r) is given. With values of 9, Q calculated 
with the guessed value of y , the radial distance (call r g ) is calculated by linear 
interpolation of the the 6 , (j) grid. Depending on the magnitude of r g - r, die next 
guess value of y is taken by the bi-section method. This procedure is repeated till 
r g « r. after getting r, 9 , <f> the x, y, z are obtained by back substitution. 
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Chapter 3 


Representation of Pelton Wheel bucket surface 


3.1 Introduction 

The primary requirement of our model is to reflect the smallest changes in the 
surface of the bucket, for which the first step is to represent the bucket surface 
with very good accuracy. 

The surface geometry of a Pelton Wheel bucket is quite complex and obvi- 
ously can not be represented by a simple analytical function. To overcome this 
difficulty the bucket surface is represented in the form 

r = f(0, (j>) (3.1) 

Here data are supplied in the form of a fine rectangular grid (M x N) in 6 and 
cf>, and at each grid point the radial distance (r) is given. To calculate the value 
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of r at any 9 and <j> between the grid points we use bi-linear interpolation. This 
interpolation is finally done (as explained in the next chapter) in the ODE solver. 
However the first problem is to create the fine grid representing r — f{9 , o) from 
a coarse grid of y — f (x, z ) data which represents the empirical measurements 
of the bucket. The fine (r, 6 , 4>) grid is created from the coarse (x, y, z) grid by 
bi-cubic spline interpolation. 

3.2 Bi-cubic spline interpolation 

We use the cubic spline interpolation scheme since it provides a superior approx- 
imation for functions that may have abrupt local changes. The objective in cubic 
splines is to derive a third order polynomial for each interval between grid points. 
In one dimension, the function would be 


fj(x) = ajX 3 4- bjX 2 + CjX + dj (3.2) 

where j(= 1,2,3 ■ ■■ ,n) 

To develop cubic splines we have to solve 4n equations for the values of ay, bj, 
Cj , dj for n intervals. These equations are obtained by imposing the conditions that 
the function yield the tabulated values at the grid points and also the piece-wise 
functions have continuous first and second derivatives at the grid points. 

The second derivatives of die above function can be represented using a first- 
order Lagrange interpolating polynomial as 
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(3.3) 


f'Ux) = f"(x j . l ) X 


Xj - 1 




■Xj 


Xj 


Xj - 1 


where /” (x) is the value of the second derivative at any point x within the j th 
(xj-i < x < Xj) interval. Integrating the equation (3.3) and using function 
equality conditions, i.e., that f{xj) must equal the tabulated fj and /(x ; _i) must 
equal fj-i, and invoking the conditions that the first derivatives at grid points must 


be equal, i.e., fj(xj) = , we get 


X J x j-i f " /„ \ , x j+i 1 c" (- \ , x j±i x j f " 


-f (xj-i) + 


3 / fo) + («j+i) 


f(xj+l) - f{Xj) f(xj) - /(Xj-i) 


^i+i x j 


Xj Xj—i 


(3.4) 


These n — 2 linear equations have n unknowns, all the f"(xj)'s. So setting 
/”(x i) and f"(x n ) to zero (which gives the so-called natural spline which has 
zero second derivatives its boundaries), the above equations can easily be solved 
by using a tridiagonal matrix algorithm. 

In two dimensional interpolation, we seek to create a fine two dimensional 
rectangular grid (M x N) of tabulated values from a smaller table of the coarse 
data To do this, we perform M one dimensional splines evaluation across the 
rows of the table, followed by one additional one-dimensional splines down each 
newly created column. This is called bi-cubic spline interpolation. 
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3.3 BHEL buckets 


Two data sets (say 1 and 2, each representing the surface of a different Pelton 
Wheel bucket) have been provided by BHEL Bhopal as test cases. We first rep- 
resent the bucket surface obtained from these data sets by the above mentioned 
interpolation scheme. 

The bucket data are provided in the form of a very coarse rectangular grid in 
x — z coordinates and y ~ f(x,z) is provided for (x, z ) grid points, along with 
the coordinates of points on the boundary of the bucket. Referring to figure (3.1), 
which shows the plan view of the buckets provided, at the points la, 2a, etc which 
lie outside the bucket, the y coordinate value does not have any physical meaning 
and so are not defined and are kept as blank entries in the input data provided by 
BHEL. This means that the y coordinate value is not defined at every' grid point 
of the rectangular grid in x — z. Implementation of bi-cubic spline interpolation 
will require complete data on a regular grid. Therefore the first task is to fill these 
blank entries of input data so that interpolation can be done. The practice used to 
fill these blank entries is as follows. 

First, the y values for the blank entries (grid points la, 2a, ..etc) are extrapolated 
by any interpolation scheme (polynomial, rational function or cubic spline) along 
the vertical z direction. If extrapolation is done across a large distance then the 
extrapolated values given by all schemes are often very unrealistic and make no 
practical sense; in such cases we abort those values and give some arbitrary' value 
(usually the y value of the neighboring grid point). This will not incur much error 
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Figure 3.1: X — z projection of the sample case of half of the Pelton Wheel bucket. 
The data provide y values at each grid point inside the demarcated boundary line. 
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in the final fine grid bucket representation, as is explained below. 

Figure (3.2) shows a bucket which is part of a sine wave ‘a’. We can see that if 
the grid point values inside the bucket are represented correctly, the interpolated 
values within the bucket are interpolated accurately (V) even if the grid points 
outside the bucket are abruptly truncated (’6’)- Thus the error in the exterior grid 
points does not affect the interpolation in the interior of the bucket, so the arbitrary 
extrapolation used does not affect the final fine-grid surface. This feature follows 
from the local nature of the bi-cubic spline interpolation scheme. The final fine 
grid surface is used for the particle tracking algori thm . 



Figure 3.2: property of spline interpolation, 
a: the actual sine curve 

bxoarse data provided for a sine curve truncated at top 

c: interpolated from b using cubic spline interpolation, showing the error intro- 
duced at outer edge does not propagate inside. 

For tire particular cases of the two BHEL data sets, after interpolation of the 
raw data in vertical direction, the data set 1 w'as interpolated from a coarse grid of 


28 




11 x 10 to a fine grid of 303 x 138, while the data set 2 was interpolated from a 
coarse grid of 10 x 8 to a fine grid of 278 x 147. The coarse grids are in x and z, 
but the fine grids are in 6 and 4> co-ordinates. This change is made for convenience 
of use in the ODE solver. The maximum limit of 6 and <f> for fine grid is calculated 
from the boundary as well as the coarse grid. 

60 

40 

20 

0 

-20 

-40 

-60 

Figure 3.3: x — z projection of the Bucket obtained after bi-cubic spline interpo- 
lation for Data Set #1 , provided by BHEL Bhopal as input data. 

Fig 3.3 and 3.4 shows the xz projection of the buckets which is obtained after 
bi-cubic spline interpolation, along with the boundary of the respective buckets. 
From the two figures (3.3 and 3.4) it is clear that the grid obtained after interpo- 
lation contains the actual bucket within it but also has points outside it. Thus it 
is necessary to identify and mark each (0, <p) grid point as being either within the 
boundary or outside it. This is done in the following manner. 

Since the fine grid is rectangular in 8 and <j>, and at every grid point r, 6 and o 
are known so the x and z coordinate for any grid point are obtained by 
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Figure 3.4: x — z projection of the bucket obtained after bi-cubic spline interpo- 
lation for Data Set # 2 , provided by BHEL Bhopal. 


x = r sin 8 cos <j> 


y = r sin 9 sin <j> 


z = r cos 6 

For each grid point, we can thus find the x and z values from (0, (p) values. For 
each x, we than find z max and z m i n , the location of the upper and lower boundaries 
by linear interpolation of the input boundary' data. If the z value of a point in the 
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bucket surface satisfies z min <z<z max , then it lies within the bucket. For each 
grid point on the fine grid we first determine whether it is within or outside the 
bucket, or on the boundary. 

We define a cell as a region which is formed by four neighboring grid point, 
i.e, cell(i, j) is formed by four grid points (i, j), (i 4- 1 , j), (z, j + 1), (i + 1, j + 1). 
If all the four grid points forming a cell lie within the bucket, then tire cell has a 
status 1, which means the cell is within the bucket. The status of cell (i,j) is 0 
if all the four grid points forming the cell are outside the boundary, so the cell is 
completely out of the bucket. The status of cell(i, j) = 2, if one or more, but not 
all four grid points, are outside the bucket; this is a boundary cell which contains 
the boundary of the bucket. 

In the particle computations, the bucket is assumed to be represented by all 
cells which have a non-zero status (if the grid is fine this will not incur much 
error). Figures ( 3.5 and 3.6) show the buckets represented by only non-zero cells. 
Thus, with this method the bucket boundary can be closely approximate to the 
exact one. The grid obtained after this procedure is now used as input to the ODE 
solver which tracks the particle motion on the bucket surface. 
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-10 0 10 20 30 40 50 60 70 80 

Figure 3.5: Bucket Grid and for Data Set # 1 along with the outer boundary pro- 
vided by BHEL Bhopal. 
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Chapter 4 


Numerical Methodology 

4.1 Introduction 

The bucket of a Pelton Wheel has been represented in the form of a fine grid in 6 
and <f>. In the governing equations (2.11 and 2.12), which track the motion of the 
particle on the bucket surface, surface derivatives at every grid point is required. 
This chapter presents the calculation of the surface derivatives and, if more than 
one bucket share the jet simultaneously, the method to find which bucket the par- 
ticle will strike. The algorithm for particle tracking, force & moment calculations 
and for the integration of many particles information to obtain the overall perfor- 
mance parameters is discussed. 
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4.2 Calculation of surface derivatives 


Our primitive equations (2.11) and (2.12) governing the particle motion on the 
bucket surface contain the terms g, §, ff, 0, gf, 0 and g§. As discussed 
earlier we have reduced the surface presentation of bucket from ip(r, 6,(p) — 0 ( 
the equation defining the bucket surface) to a r - f(9,<j>) =0 form. So g = 


X > dr* 


n AAL — n SA — _££. §± — _<LL <¥± _ 3 2 / 

u ’ grd0 u ’ 89 86' dip 8<t> ’ <50 2 ~ ’ 


d<p 2 


-0 and 


- a 2 / 

868 $ 868 $' 

The value of r as a function f(9, (j>) of 6 and cf> is represented by the fine grid 
generated by the bi-cubic spline interpolation procedure. We now need to use 
this data to calculate the derivative of / at any arbitrary (#, <j>) point which the 
particle may be at during its journey on the bucket surface. To ensure continuity 
of derivatives of the interpolated function across the cell boundaries, the value of 
first-order, and second-order derivatives with respect to 9 and <f), is obtained at the 
grid points using a central difference scheme and the value of the function and 
these derivatives at any arbitrary point is obtained by linearly interpolating each 
of these separately. This procedure is elucidated below. 

The function value ( / (6, 4>) = r ) is available at every grid point on the fine 
grid. We compute the first and second order derivatives at the grid points using 
central difference, i.e. 


d£ 

d9 


8i,$j — fd(9i, (pj) — 


2A9 


and similarly, 


(4.1) 
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— le , — f (Q. a.) — IMtA 
dp 9 "** 2A , 


(4.2) 


Similarly the second order derivatives at the grid points are computed as fol- 
lows: 


av 

de 2 



= fee(9i) 4>j) 


/ <j>j) ~ 2/ (flj, 0j) + / (0j-i, (ftj) 
A0 2 


(43) 


and similarly, 


a 2 / 

df 


6i,4>j ~ f<P<i>{Qii 4>j) — 


f{Qu ‘ftj+i) ~ 2/ (6>i, </>■,-) + f(Qu ) 


(4.4) 


and the cross-derivatives at the grid-points are computed by 


d 2 f , m+b ^--n) + m-t> fr-o - m+u ^,-i) - /(^-ir <>j-n ) 

ddd<l> ^ ~ i9 *^ u ~ 4A0A?i 

(4.5) 

At the boundary points a forward difference or a backward difference is used in place of die cen- 
tral difference scheme. 

To find the value of the function and its derivatives at the co-ordinates (0. <j>), 
we first find its cell which is bounded by the grid points (0*, 4>j), (0i+i , 4>j), {9u A>+ 1 ) 
and (0 i+1 , <l>j+i), then linear interpolation is used 
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/(M)= m,«(%^)(^) + 

/(«<+!■« (^) + 
m. *j+i) (^Sr) ( ^ ■■ 1 *y-»> ) + 

/(»«• i,*+i) ( ^fc^a ) ( ^fcfeirfi ) 

(4-6) 


where / could be / , f e , f^, f gs , f^, This interpolation method is fast and re- 

liable and ensures continuous derivatives across cell boundaries. This last fac- 
tor is very important, as the particle dynamics can be in error if the derivatives are dis- 
continuous. 


4.3 Algorithm for multiple buckets 

If the number of buckets in the proposed Pelton wheel design is not already known 
then the analyses can be carried out for single bucket. However the program 
is generalized to handle the multi-bucket case, where we know the number of 
buckets involved. 

Figure (4.1) shows the schematic diagram of a jet with buckets as in a practical 
Pelton Wheel system. The jet is shared by three buckets at a time. In this figure, 
The operating range of angles for the given geometry is Xmin to Xmai which 
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Figure 4. 1 : Schematic diagram of jet along with the buckets to which it is striking. 
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are 


Xmin — ~ arccos ((/l — l 0 ) J Ri ) 


and 


Xmax = arccos ((h - Iq)/Ri) 

for a single bucket, where Xmin and Xmax are the angles at which contact between 
the jet and the bucket is started and finished respectively. 

When a number of buckets are sharing the jet simultaneously, then for pth 
bucket 


Xpmin — Xmin (p 1) X 2tt/(J 

where p=l,2, or 3 etc., indic|te the first, second and third buckets in the ‘line of 
fine’ of the jet. Similarly 


Xpmax — Tflin ( Xpmin "h 271"/ C[, Xmax) 

where q is total number of buckets in the Pelton Wheel system. The pth bucket 
starts its motion from Xpmin at t = 0 reaching xpmax at t = n x dt, where n is 
no of particles [n = ( Xpmax ~ Xpmin) /(w * dt)} and dt is the time step. Starting 
from t = 0, particles are assumed to come at regular intervals of dt from random 
positions in the cross-section of jet. 
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If a single bucket is moving, then within its range of operating angle, the whole 
of jet may not strike the bucket and part (depending on the geometry of the bucket) 
of it may go to waste. As the particles are coming out from random positions in 
the jet cross-section, therefore it is necessary to know whether the particle co min g 
out at any instant of time strikes the bucket. The contribution of all particles which 
do not strike the bucket is made zero and their further investigation is not done. 
Generally the number of buckets in a Pelton Wheel turbine are of the order of 
15 -20, therefore at a time the jet is shared by two or more buckets (seldom > 
3). Each particle coming from the jet may strike any one of the buckets sharing 
the jet at that particular instant of time. Since the angular position of buckets are 
different, the force contribution by the particle will depend on which bucket it 
strikes. Therefore it is necessary to know which bucket, it strikes. This is checked 
as follows. 

The angle of the arm with the Z axis when pth particle strike at time t v is 
X P , the coordinates of the particle with respect to center line of the jet are AX P 
and A Z p . From figures (4.1 and 2.1), the projection of qth bucket on XZ plane 
at time of impact of pth particle = (Ri — R^) cos x Pq cos a, where a is the angle 
of inclination of bucket with its local frame of reference. The Z coordinate of 
particle with respect to the fixed frame is h+AZ p 

If 


R 2 cos Xpi < h + A Z p < R\ cos x P i (4-7) 
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is not satisfied then the particle will not strike the first bucket, 
or if 


Ri cos x P i <h + A Z p < R x cos x v i (4.8) 

is not satisfied along with equation (4.7) then the particle will not strike the second 
bucket also and will strike to third bucket (note that here xp 2 =Xpi+27r/g etc. The 
idea can be extended to more buckets if more share the jet at any time). 

However the above criteria considers the vertical ( Z coordinate) direction 
only, and only assures that the particle will not strike the bucket if the criteria are 
not satisfied. Fulfilling the above criteria does not assure that particle will strike 
the bucket as it may pass from the side of the bucket. If the above condition is 
satisfied for the particle striking on pth bucket with coordinates X p and Z p , solve 
the equation 2.33 by the scheme mentioned in the second chapter. If X p and Z v 
are such that they are out of the projection of the bucket, the equation will not 
converge. If the equation does not give a converged solution in a limited number 
of iterations (say 1000) then the particle is assumed to have missed the bucket and 
is arbitrarily assigned values of 9 and tp which must be greater than 9 rnax and (p max . 
As already mention in chapter 3, the bucket is represented by only the cells which 
have non-zero status. If 9 min < 9 < 9 max and (p m i n < (f> < rpmax, and the 
cell to which 9 and <p belongs have a non-zero status, then the particle is assumed 
to strike the bucket, otherwise not. If the particle does not strike the bucket then 
its efficiency of conversion is made zero and its path is not computed. 
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4.3.1 RK4 solver 


The equations (2. 1 1) and (2. 12) track the motion of the particle over the bucket 
surface, and are solved as an initial value problem requiring four initial condition 
in terms of 9, dp, 6,<j> which are obtained as explained in the second chapter. With 
these initial conditions the value of 9,<j>, 9,<f> at subsequent time steps is calculated 
by a fourth order Runge-Kutta ODE solver. After each time step we check whether 
the particle is within bucket or not. The particle can leave the bucket in one of 
these two ways: 

1. If the status of the cell with the new time step values of 9 and <j> is zero (i.e 
it is out of bucket), or 

2. if the normal acceleration (and hence the force of the particle with respect 
to bucket) is negative, (i.e., bucket pulls the particle); this means that the 
particle would lose contact with the bucket. 

In applying the second criteria we assume that the particles are non-interacting, 
But in reality this is not true, because the flow is continuous; so if any single 
particle’s acceleration becomes negative it will try to leave the bucket but will be 
pushed back by the main stream around it, at the expense of the kinetic energy of 
the main stream. So two different sets of criteria are applied to check the presence 
of the particle at each time step 

case 1 : The particles are considered as non-interacting and both conditions 
(1,2) mentioned above are checked. If either occurs, the particle leaves the bucket. 
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case 2 : Every particle which strikes remains in the bucket till the boundary, 
thus only tire first (1) condition is checked. 

4.3.2 Calculation of performance parameters 

As already discussed, the ODE solver gives the value of 64 , 64 at the next time 
step. The normal acceleration of the particle is given by 



The velocity and acceleration are computed by calculating the change in the 
94 , 64 with time, the kinetic energy of each particle is calculated at every time 
step. The difference between the initial kinetic energy and kinetic energy at exit 
is used to calculate power (the energy transferred to the bucket) based on kinetic 
energy approach. It is assumed that there is no energy loss during the motion on 
the bucket surface, as it is frictionless. In the torque approach the torque (2.25 and 
2.26) produced by each particle present in the bucket and the cumulative torque at 
every time step is calculated after calculating the acceleration of the particle at that 
time step. The cumulative torque (T =E t r , where r is the total torque produced 
at any time step by all particles present in the bucket at that time step) produced 
for a given bucket angular velocity is used to calculated the energy transmitted. 

As already discussed, the particles are represented by the random numbers in 
the jet cross-section. Therefore, to minimize the effect of randomness' in the re- 
sults, the computations are done for large number of separate simulations. Each 
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simulation has a separate set of random numbers and hence a separate set of initial 
positions on the bucket surface where the particles strike. For each simulation, the 
performance is calculated in the way mentioned above, and finally all the perfor- 
mance parameters are calculated by averaging the results for different simulations. 
The use of different criteria (for loss of contact) yield slightly different values for 
the bucket performance parameters, as will be seen in the next chapter. 
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Chapter 5 


Results and Discussion 


In tills chapter, the performance parameters are calculated for a test case for which 
the exact analytical solutions are known and thereafter the simulations are done 
for the two real Pelton Wheel buckets provided by BHEL Bhopal. 

5.1 3 D Plane Bucket 

The surface over which we study the motion of the particle is a planer surface 
ABCD as shown in figure (5.1), this surface rotates around the X-axis. The co- 
ordinate system ( xyz ) rotates with the plane, and the co-ordinates system (A YZ) 
is fixed in the lab frame. We also use the spherical co-ordinates (r,0, o) in the 
rotating frame with origin o'. The value of r at any arbitrary point ( 6 , <j > ) is gener- 
ated using tlie interpolation scheme. The jet strikes the plane at a height h, coming 
parallel to Y axis, the jet transmit the power to the plane due to the initial impact, 
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but the particles then lose contact with the plane. The analytical solution was done 
for the test case by Gaurav Dar[7]. The comparison of the results by the program 
with analytical results are shown in table (5.1). it is clear that the test case has 
good match with exact solution. 



Figure 5.1: Planer bucket in 3D 


5.1.1 Analytical solution 

The equation for the surface in fig (5. 1) with respect to the co-ordinate system 
{x, y, z) can be written as y = z + 1, here tire value of R is taken equal to 1, tliis 
equation can also be written in the fonn 
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ip = r : — — , 

sm 9 sm <p — cos 9 

P is the position of the particle on the bucket surface, its distance from X axis 
is p, which is given by 

, o o, 1 r , o ~ . /•> 9r\ 2 COS 9 1,1 

p = (y 2 + z ) 2 = r[sin^ 6 sm ^ + cos 2 9 -\ — — — h — ] 3 


The total torque generated by the particle during its motion starting from t = 0 to 
t = tf is 


T = — 2 J 5 (coppX + 2 tupXpJ dt 


(5.1) 


and the total torque exerted along the X-axis during the same time interval is 


Tx = -2 / uppdt (5-2) 

J o 

The power delivered by the particle on the bucket is given by 

P = - J T x udt ( 5 - 3 ) 

P and A KE (change in kinetic energy) should be equal. 
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Table 5.1: Test case (3 D plane bucket) for the validation of code. 


Run 

h 

a rs-impact 
exact 

pimpact 
r exact 

A K im P act 

y num 

pimpact 

1 num 

T1 

1.0 

o.2i x io- 3 

0.12xl0 -3 

0.21 xlO" 3 

0.12xl0- 3 

T2 

2.0 

0.09 xl0~ 3 

0.08 xl0“ 3 

0.09 xlO" 3 

~0.08xl0 -3 


The numerical results are compared with the analytical results for the param- 
eters w,v, h, where w (angular velocity of bucket) = 0.2, v (vel of jet) =1.0, h 
is the height of jet above the .pivot, P = power due to impact ; A K- kinetic 
energy lost at impact and P in = mv]dt / 2 (initial kinetic energy of particle), 
^j^impact _ p K a fter impact- The nutial conditions (^Oj^Oj^o) are 7r/2,7r/2 
and -0.5 respectively for both test cases, and </> 0 is 0. and 0.5 for test case 1 aid 2 
respectively. 


5.2 BHEL buckets 

In the following simulations the performance parameters (efficiency, torque and 
the power output) for the real buckets are found out for 20 separate simulations 
and then averaged. While calculating for performance parameters both criteria 
(for loss of contact) as mentioned in the previous chapter are used. Two separate 
simulations are carried out, one for the single bucket and another for the multi- 
bucket case where the number of buckets in the system is given as 20. The results 
obtained for both cases are tabulated below. 
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5.2.1 Single bucket case 

Only one bucket is present in the line of action of jet, the simulation is carried out 
for the whole span of the bucket in which the jet can strike it. 

For Data Set # 1 

time step (dt) = 0.0001 sec 

Xmin =-0.7935, Xmax =0.7935 (entire range of striking angle) 

#1=222. 12 mm, #2=126.92 mm (figure 4.1) 
angular velocity (u) =86.5 rad/sec 
velocity of jet (u,)=31 M/sec 
radius of jet (i 0 )=20.25 mm 

inclination angle (a) = 8° (rotation of bucket about center of local frame of 
reference) 

total number of particles (n) = 184 
For Data Set # 2 
time step (dt) = 0.0001 sec 
Xmin =-0.76178, Xmax =0.76178 
#1 = 216.41 mm, # 2 =124.58 mm 
inclination angle (a) =4.77° 
total number of particles(n) =178 

other input parameters for data set # 2 are same as for data set #1 

criteria 1 : Particles are non interacting. 

criteria 2 : Particles can leave the bucket only at boundary. 

Table (5.2) shows the results for the two buckets. It is clear that first bucket 
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Table 5.2: The average results for the 20 simulations, for the two loss of contact 
criteria ( 1 for non interacting particles, and 2 for particles leaves only at bound- 
ary). Here rji is the efficiency based on how many particles actually strike the 
bucket and % is the overall efficiencies based on the total number of particles. 




Criteria 

Vi(%) 


1 

1 

1 

78.52 

60.75 

2 

2 

1 

59.58 

44.18 

3 

1 

2 

79.31 

62.01 

4 

2 

2 

72.77 

53.96 


is more efficient than the second one from both criteria. Recall that positions of 
the fluid particles in the jet is random and there is only one bucket, hence not all 
fluid particles need to actually strike the bucket. Note that in the single bucket 
case many particles miss the bucket and r/ 2 is substantially lower than rji (which 
is obtained only for particles which are striking). 


5.2.2 Multiple buckets 

In a real life Pelton Wheel system, the jet is always shared by two or more buckets 
at any time (seldom > 3). We have simulated 3 buckets simultaneously sharing 
the jet (total number of bucket =20). Simulations are done for each bucket indi- 
vidually and the overall performance is calculated by summing up the individual 
contributions. In chapter 2. it is mentioned that energy transfer by the particle is 
of two kinds: first, during impact there is loss of kinetic energy and part of the 
kinetic energy is converted into power (call it impact power), and secondly power 
transmission occurs as the particles move on the surface of the bucket (i.e power 
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due to motion). In addition to the overall efficiency, we also compute two other 
efficiencies rj impact and ri motion . The first of these tells about the loss (averaged 
over all particles) of the initial kinetic energy due to the particles’ impact with the 
buckets. Thus 

’ K.E after impact +] 30 wer transmitted to bucket on impact 

/impact Initial Kinetic Energy 

The r) mot i on tells about the amount of kinetic energy lost because of the kinetic 
energy carried out by the particles finally exiting the bucket. 

K.E after impact— exit K.E loss 

i /motion — K.E. after impact 

These two efficiencies give different perspectives on the bucket performance: 
Vimpact is used to check whether the operating conditions are optimum or not and 
Vmotion basically tells about the optimization of the surface profile. Tables (5.3 to 
5.6) shows the results for the two different buckets with different criteria of loss 
of contact. It is clear from there that rj motion is generally lower showing that exit 
loss are significantly decreasing bucket performance. 

To calculate power due to motion we use both the kinetic energy and torque 
approaches. In the first method, we find power transmitted during motion from the 
change in kinetic energy of the particles (as motion is assumed to be frictionless). 
In the second approach, the torque applied by the particles is computed and the 
power is computed from the torque and the angular speed. The two different 
methods offer a vital check on the accuracy of the computations, as they should 
give identical results. It will be shown further that both approaches yield almost 
same results in our computations, both in terms of numeric values as well as the 
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qualitative trends. The results for multi-bucket case for both data sets provided 
by BHEL are tabulated below. In all the cases hereafter, the total efficiencies are 
based only on the striking particles (as almost all particles strike some bucket). 
For both data sets 
time step(df) =0.00003 sec 
total number of particles(n)=121 

total number of buckets =20, rest of the input parameters are same as in their 
respective previous cases for single bucket analyses. 


Table 5.3: Data Set 1 Criteria 1 


Bucket number 

1 

2 

3 

Total 

r?(%) 

47.84 

40.38 

0.07 

88.29 

V*(%) 

46.13 

38.8 

0.07 

85.0 

Vimpacti %) 

53.10 

43.41 

0.08 

96.59 

^motion (%) 

44.31 

36.15 

0.065 

80.53 

Impact Power 

0.807 



1.848 

Impact KE loss 

1.0362 

1.4645 

0.0022 

2.503 

Power! due to motion 

8.732 

6.7083 

0.013 

15.453 

Power 2 due to motion 

8.044 

6.4069 

0.013 

14.464 

Avg no of particles strike 

65.7 

55.2 

0.1 

121 


In tables (5.3 to 5.6) rj and rj* are the efficiencies calculated based on kinetic 
energy and torque approaches respectively, subscript 1 is used for power calcu- 
lated from kinetic energy, and 2 is for torque. All energy and powers are in kW. 

The tables (5.3 to 5.6) shows that there is generally a close match between die 
parameters calculated based on die kinetic energy and die torque approaches. This 
means the computations are quite accurate (except perhaps in Table 5.4). Also 
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Table 5.4: Data Set 2, Criteria 1 


Bucket number 

1 

2 

3 

Total 

t?(%) 

38.75 

38.15 

0.38 

77.28 

??*(%) 

34.82 

35.66 

0.38 

70.86 

r ?impact( % ) 

49.25 

43.66 

0.66 

93.57 


30.71 

30.19 

0.2 

61.1 

Impact Power 

1.3324 



3.22 

Impact KE loss 

1.769 




Power! due to motion 

6.104 

EE2I 

0.038 

11.61 

Power 2 due to motion 

5.349 



10.732 

Avg no of particles strike 

62.35 

57.75 

0.9 

121 


Table 5.5: Data Set 1, Criteria 2 


Bucket number 

1 

2 

3 

Total 

77(%) 


41.46 

0.07 

92.65 


51.38 

41.62 


ESEBI 



43.41 

0.08 

96.59 

^motion.^ 

49.87 

39.2 

0.065 

89.13 

Impact Power 

0.807 




Impact KE loss 

1.0362 

1.4645 



Poweri due to motion 

9.003 

6.916 



Power 2 due to motion 

9.0528 

6.947 

0.013 

mm 

Avg no of particles strike 

65.7 

55.2 

0.1 

121 
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Table 5.6: Data Set2, Criteria 2 


Bucket number 

1 



2 

3 

Total 

v(%) 

44.92 

41.26 


86.77 

??*(%) 



0.59 

86.01 

^impact^) 




93.42 


41.03 

36.65 

0.41 

78.09 

Impact Power 



0.035 

3.22 

Impact KE loss 



0.051 

4.482 

Powerx due to motion 


mi 

0.078 


Power 2 due to motion 

7.148 

6.056 

0.078 


Avg no of particles strike 

62.35 

57.75 

0.9 

121 


note that criteria 2 (which seems to be more practical) gives higher efficiencies 
for the same configuration. Also, the power due to impact is very low and most 
of the energy conversion takes place during motion of the particle over the bucket 
surface. 

A perfect bucket profile (ideal case) is one in which there is no energy loss at 
impact (i.e the difference between jet velocity normal to bucket and velocity of 
bucket in normal direction at the point of impact is zero) then the only energy loss 
will be the energy carried away by the particle at exit, but real cases are always 
associated with impact losses and thus lowering the overall efficiency. Therefore 
efforts are made to reduce the kinetic energy loss at impact by changing the con- 
figuration of the system (like height of jet above pivot, pitch circle radius, angular 
velocity of bucket etc.) with an objective to get the optimum operating conditions 
for a required power output under given head. Similarly the kinetic energy loss at 
exit can be reduced by changing the bucket profile. With this program the effect 
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of changing various parameters as well as the bucket profile on impact kinetic 
energy loss and the power due to motion can be studied to find out the optimum 
operating conditions and the bucket surface profile. 

The selection of time step (dt ) of advancement for ODE solver is also impor- 
tant. As the number of particles depend on time step; the smaller the time step, the 
larger will be the number of particles being used to simulate the flow. Table (5.7) 
shows the efficiency variation with time step size for data set # 1 with criteria 2. 
The results are almost same demonstrating that our results are robust. 


Table 5.7: Efficiency for Data Set #1 for different time steps of advancement for 
the ODE solver 


dt / 7](% ) 

1st Bucket 

2nd Bucket 

3rd Bucket 

Total 

0.0001 

52.71 

40.26 

- 

92.97 

0.00003 

51.38 

41.46 

0.07 

93.07 

0.00001 

52.03 

40.98 

0.08 

93.09 

0.000005 

54.58 

38.35 

0.18 

93.11 


Since the jet is modeled as a series of non-interacting plastic particles, each 
particle has its own independent time history and hence, the efficiency of conver- 
sion for each particle is different and depends on the angle of arm with vertical at 
the time of impact and its position on bucket surface where it makes the impact. 
Figure (5.2) shows the variation of the efficiency of the particles with respect to 
the angle of bucket arm with vertical when the particle strikes. Here only one 
bucket is moving through its full range of operation. It is clear from this figure 
that the efficiency generally decreases as the rotational angle increases, and the 
particles which are striking in the last phase of the contact have extremely low 
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efficiencies thus lowering the average efficiency. Therefore it is not advisable to 
let the bucket operate in its full range. Clearly the optimum range of average ef- 
ficiency can be selected from the figure, which will then give the range of angle 
of operation for each bucket and the number of buckets which will share the jet 
at any time. Hence the total number of buckets required in the Pelton Wheel can 
be chosen after fixing the efficiency required. Also larger the number of bucket, 
the higher will be the average efficiency; but cost will increase, so a compromise 
between efficiency and the cost can be made. 



Figure 5.2: Variation of efficiency with rotational angle (jet ratio=8.75) 

1. data set #1 

2. data set #2 

certain 

The jet ratio (ratio of pitch circle diameter of Wheel to jet diameter) in the 
Pelton wheel varies from 5 to 20, the efficiency of the Pelton Wheel turbine is 
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also depencertaindent on the jet ratio. Figures (5.3 to 5.6) shows the efficiency 
variation with rotational angle with different jet ratio. It is obvious from the fig- 
ures that for higher jet ratio (i.e low jet radius), not only is the average efficiency 
somewhat higher but the curve for efficiency has a lower slope up to a certain ro- 
tational angle. This effect is however less prominent beyond a certain jet ratio (in 
this above 8.75). As the higher jet ratio means lower flow rate and hence less total 
power output but with higher jet ratio, the maximum number of jets (which can 
be used with a particular configuration) and hence the power output will increase, 
therefore the jet ratio should be optimised from efficiency, power output and cost 
considerations. The present simulations are done (wherever not mentioned) for jet 
ratio of 8.75. 



Figure 5.3: Efficiency vs rotational angle (jet ratio =17.7) 

If the number of buckets for the proposed design whose performance is to be 
computed is fixed earlier (from past experience or due to any other constraints) 
then from the geometry we can know what is the range of operation of a single 
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Figure 5.4: Efficiency vs rotational angle (jet ratio= 11.82) 



Figure 5.5: Efficiency vs rotational (jet ratio=5.9) 
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Figure 5.6: Efficiency vs rotational angle (jet ratio=5) 

at any time. For a variety of designs it is seen that the number of buckets that 
share the jet at any time is rarely more than three, this computer model at present 
computes the case for three buckets simultaneously, but it can be easily extended 
to as many as required. Figures (5.7) and (5.8) shows the efficiency variation for 
tire two buckets whose results have been discussed, here three buckets are sharing 
the jet simultaneously, but as clear from the tables that the contribution of the third 
bucket is very low in comparison with the other two buckets, note that the average 
efficiency is increased up to a large extent. 

When jet is simultaneously striking more than one bucket, each bucket will 
produce a different torque that depends on the time history of the particle. The 
total torque and hence the power will be calculated by su mmin g up the individual 
contribution of each bucket. Figure (5.9 and 5.10) shows the variation of torque 
widi time for each bucket that is sharing the jet for both data set #1 and data set #2 
respectively. The total torque variation for both buckets is shown in figure 5.11. 
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0.8 


0.75 » 1 1 1 1 1 1 1 

-0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 

rotational angle 

Figure 5.7: Efficiency vs rotational angle for Data Set#l 
a. based on torque; b. based on kinetic energy 



-0.8 -0.7 -0.S -0.5 -0.4 -0.3 -0.2 -0.1 

, rotational angle 

Figure 5.8: Efficiency vs rotational angle for Data Set 2 
a. based on torque; b. based on kinetic energy 
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This figure shows that the bucket is free of fluid within 0.01 seconds. With this 
knowledge, we can decide how many jets can be used in a multi -jet Pel ton Wheel. 



Figure 5.9: Variation of total torque with time for data set #1 
1 , 2 are the the buckets which are sharing the jet 

It has been already discussed that each particle has its own independent time 
history. Thus the path of each particle can be traced while moving over the bucket 
surface. Figures (5.12) and (5.13) show the path of few particles with their indi- 
vidual efficiency of conversion (criteria 2) for data set 1 and 2 respectively. Since 
the efficiency of conversion of each particle is known, by tracking the path of 
particles over the bucket surface we can locate the regions in the bucket surface 
which give low efficiencies. 

The path of particles over the bucket surface using criteria 1 (non interacting 
particles) shows (figures 5.14 and 5.15) that some particles leave the bucket in the 
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torque 


Figure 5.10: Torque vs time for data set #2 
1, 2, 3 are the three buckets to which jet is striking simultaneously. 



Figure 5.11: Total torque vs time * 


1. DataSet# 1 

2. Data Set # 2 





middle. This indicates local irregularities in the surface of bucket at those places 
where the particle has left the bucket, by modifying these areas the bucket surface 
can be optimized. For a perfect bucket surface the performance of the bucket 
with both criteria should be equal, so closer the performances of the bucket with 
two different criteria mentioned above, nearer the bucket surface will be to the 
optimum one. 



•10 0 10 20 30 40 50 60 70 80 

xaxis 


Figure 5.12: Path of few particles with their efficiency while moving over the 
bucket surface obtained for data set #1, criteria 2 (jet ratio=5) 


The results for the above mentioned buckets have been found to have close 
match with the experimental results of BHEL Bhopal. The program was tested by 
varying the input parameters, the trend as well as the prediction of the performance 
parameters were found to be on the right track. 
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xaxis 

Figure 5.15: Path of few particles with their efficiency while moving over the 
bucket surface obtained for data set #2, criteria 1 (jet ratio=5) 

5.3 Conclusion and future work 

In this thesis, the program which started by Ajay Kumar and followed by Gaurav 
Dar, is now fully developed. One test case (also tested by Gaurav Dar during 
development stage of the program) is simulated with the completely developed 
program, and finally the performance parameters for two Pelton Wheel buckets 
are computed and compared. The complete program is developed to handle a 
variety of cases and to simulate the effect of any of the input parameters, the 
present program is tested against large set of experimental results and found to 
be quite sensitive to the surface changes which was the prime requirement. Thus 
this program is extremely useful for industrial application. It can act as a tool to a 
designer by which he can calculate the performance parameters of a Pelton Wheel 
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numerically These parameters includes 


1. The overall efficiency of the bucket. 

2. The kinetic energy loss associated with impact ( rjimpact )• 

3. The kinetic energy exit loss ( 7] m0t i 0 n )• 

4. The efficiency vs arm angle (x) curve, which tells at what striking angles 
the performance is good. 

5. The particle paths can be used to correlate lower efficiency obtained (for 
certain angles) with portions of the bucket surface. This will help in locating 
the less optimum portions of the bucket surface. 

6. Total power output of the complete Pelton Wheel. 

7. Torque variation of individual buckets and the total torque variation with 
time. 

8. Number of jets which can be used with given configuration and bucket sur- 
face profile. 

The program gives the liberty to the designer to think of many different designs for 
a particular requirement, since different designs can be tested and compared with 
the help of the program and the most efficient can be found. The most efficient 
design can be further optimized by analyses of particle’s individual efficiency and 
its path over the bucket surface. Also the most optimum operating conditions for 
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the required output can be found out by changing the operating conditions and 
studying its effect on the amount of kinetic energy loss at impact. 

In the future the present work can be extended for non-friction surface. The 
further work can include the vibration of rotor and its effect on perfonnance of 
the Wheel. Attention can also be paid upon the water jet which is assumed to 
be running in straight line while in reality it has some curvature. Here the fluid 
coming out from the jet is considered as a cylindrical stream, while in practical 
cases due to the presence of a spear rod the fluid flow is like flow between two 
concentric cylinders where the inner cylinder radius is quite small. 
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